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Executive Summary 

The goal of this project was to col- 
lect “research software” written in 
IDL 1 to support analysis of IRAS 
data and make it available to the 
larger community. “Research soft- 
ware” describes software created by 
researchers and staff for a specific re- 
search goal, but lacks sufficient doc- 
umentation, easy-to-use interfaces, 
and rigorous debugging. Additional- 
ly, most of the IDL/IRAS code avail- 
able needed to be ported to a 
(largely) hardware independent new 
version of IDL. This report summa- 
rizes the activities and results of the 
one year grant. Due to the absence of 
Prof. Domik during the University 
year 1990/1991, the project received 
a no-cost extension of six months. 


1 . “Interactive Data Language’ by Re- 
search Systems, Inc. 


This report details: 

• a list of IRAS/IDL functions that 
were created/modified under this 
grant; 

• examples of these functions. 

All software modules were incorpo- 
rated into the IDL Astronomy Li- 
brary by W. Landsman at GSFC. 
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IRAS/IDL functions created/modified 


Most of the software listed below was produced over the last five years by re- 
searchers, students, staff or visitors at CASA other than proposal investigator and 
research assistant An exception is the software to preprocess IRAS images from 
skyflux data, which was developed by Gitta Domik. All other programs were mod- 
ified, documented and tested by C. Scott Merkle. 

All functions are now available through ASTRON.TLB, an astronomical IDL li- 
brary, supported through a different NASA grant and supervised by Wayne Lands- 
man at GSFC. 


Programs Description 

IRAS_DISPLAY Displays one selected band of IRAS skyflux 

plate. Several display options are available. 


Procedures Used by IRASJDISPLAY 


COLORBAR 

Displays color bar on die right edge of the 
display. 

CONT_HARD 

Creates file of contour plot of IRAS image that 
can be printed on a postscript printer. 

CONT_OVER 

Overlays a contour plot of an image on the 
image display. 

CURS_BLOW 

Zooms a section of an image. 

DISPLJSTRETCH 

Converts a floating point image to a byte image 
giving the statistically best min/max range. 

DWEK 

Finds the Eli Dwek temperature from the ratio 
of two IRAS bands. 

IRAS AXIS 

Draws the coordinate RA and DEC axes on the 
screen corresponding to the appropriate plate. 
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IRASBHR 

IRASHDR 


Obtains plate header information from the 
binary .bhr file 

Obtains plate header information from the 
ASCH .hdr file. 


iras_read 

irasjreset 
^ IRASRAXY 

IRASXYRA 

STRBLK 

STRNAM 


Obtains a floating point array of an IRAS 
image from disk. 

Resets the graphics system variables. 

Obtains the sample andjjne number for a 
given IRAS plate and specific RA and DEC. 

Obtains the RA and DEC in both decimal and 
h,m,s for a given IRAS plate no. and line and 
sample number. 

Removes all the blanks from a string. 

jt 

Removes aU the characters that cannot be used 
in a filename. 


Supplementary Programs for IRAS processing 

ir 

AP MASK Calculates the subanray used to determine 

the contribution of ea*$ pixel in a box array 
to a simulated circular aperture.(Used by 



IRAS_FLUX). 

BANDS 

Obtains floating point and calculates byte 
arrays of all the bands of one IRAS HCON 


image. 

BB_FLUX 

Gives the blackbody flux. 

CONTAXIS 

Labels RA and DEC axis on a contour plot 
drawn by CONT. 

DECIMAL 

Converts from degrees to decimal. 
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GSJMAGE 

HCONS 

IRASCENTER 

ERASLST 

IRAS_STRETCH 

IRAS.CURSOR 

IRAS_FEND 

IRAS_FLUX 

IRAS_GS_QMS 

ra_blowup 

STARJPLOT 

STRDEC 

TMPR 

QMS_CT 


Scales an image according to a color table. 

Obtains floating point and calculates byte 
arrays of all IRAS HCONS of one 
wavelength(band) . 

Calculates the center RA and DEC for a given 
IRAS plate. 

Assembles header file from FITS header 
information. 

Linear stretching/compressing of floating point values 
to byte values. 

Obtains the RA and DEC for a selected 
point on an IRAS plate. 

Obtains the plate and tape information given 
a specific RA and DEC. 

Obtains RA, DEC, and flux value for a selected 
point on an IRAS plate. 

Creates an image file to be printed on a QMS 
printer. 

Gets a section of an image using RA and DEC 
input and zooms it 

Marks the star type on a contour plot given 
the RA,DEC and Star Type. 

Converts a string array, where each element 
has 3 components (hr,min,sec), to a floating 
point array. 

Calculates the color temperature of an IRAS 
image. 

Loads a color table that has the same scaling 
as IRAS_GS_QMS . 
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WINDOW_TO_PS Reads a window from an X-window display 

saving it to an array and to a postscript 
formatted file for later printing. 

Preprocessing tools for IRAS skyflux plates 

IRAS_FLATTEN Flatten the background of an IRAS image. 

Procedures used by IRAS_Flatten 

PLANE Fits a plane through 4 points(least squares 

fits) 

HP Fits a hyperbolic paraboloid through 4 points. 

IRAS DESTRIPE Removes periodic stripes from the IRAS image. 

Procedures used by IRAS_Destripe 

FRQ Estimates the frequency of stripes in IRAS 

skyflux image. 

CALJSLOPE Estimates the slope of stripes in IRAS skyflux 

images. 

Merging and mosaicing routines 

ADDPLATE_3 Adds or merges three hcons of an IRAS plate 

ADDPLATE_2 Adds or merges two hcons of an IRAS plate 

ADD_BANDS Reads the floating point and byte arrays of all 

the bands of an addplate image. 
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PLATES 


Combines two IRAS plates that are next to each 
other and have the same declination or radial 
ascension. 


Procedures used by PLATES 


HORZ 


Combines two IRAS plates that are next to each 
other and have the same declination. 


IMGET 

MINMAX2 

VERT 


Reads the floating point array of an IRAS 
image from a .dat file. 

Converts floating point image to byte non- 
statistically. 

Combines two IRAS plates that are next to each 
other and have the same RA. 


Examples 


The next pages show samples of IDL/IRAS functions. 


Publications describing efforts of this grant: 

Domik, G., Brugel, E.W., Stencel, R.E., Vasudevan, S. t Pang, J. 1990, 
Workstation based Preprocessing of IRAS Skyflux Images, Publications of the 
Astronomical Society of the Pacific, October 1990. 

Domik, G., Brugel, E.W., Stencel, R.E., Vasudevan, S., Pang, J. 1990, Applica- 
tions of IRAS Preprocessing at the Workstation, 176th meeting of the American 
Astronomical Society, BAAS Vol. 22, No. 2, 1990 
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List of Figures 


Figure 1: Use of IRAS_DISPLAY to display Skyflux Plate 28 with coordinate 
grid. Andromeda galaxy is to the lower right 

Figure 2: IRAS_DISPLAY offers optional annotations, such as contour lines. 
Figure 3: Zooming up Andromeda galaxy. 

Figure 4: Ratio images calculating the Eli Dwek temperature can be created 
through the program DWEK. 

Figure 5 (a): The original Skyflux Plate 75 (band 3) shows sensor deficiencies in 
need of image restoration. 

Figure 5 (b): Skyflux Plate 75, band 3, after flattening (IR^_FLATTEN). 

Figure 5 (c): Skyflux Plate 75, band 3, after flattening and destriping 
(IRAS_DESTRIPE). 

Figure 6: Result of mosaicing several IRAS Skyflux Plates (PLATES): 
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In the center is Lambda Orionis. The color table from Figure 1 is used, after a la 
statistical stretch on the original skyflux values. 

Figure 7: To the left is HCON 3 (band 4) of plate 28. The missing scanlines visible 
in the Andromeda region are filled out by data from HCON 1-2 in the right image 
(ADDPLATE_3). The scale at the bottom relates the pseudo color of the image to 
low data values (blue) up to high data values (yellow). 
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Figure 1: Use of IRAS.DISPLAY to display Skyflux Plate 28 with coordinate 

grid. Andromeda galaxy is to the lower right c&iGINAL PAGE IS 

OF POO* QUALITY 





Figure 3: Zooming up Andromeda galaxy. 
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through the program DWEK. 







e 5 (a): The original Skyflux Plate 75 (band 3) shows sensor deficiencies in 
of image restoration. 






Figure 6: Result of mosaicing sever al IRAS Sky flux Plates (PLATES): 
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In the center is Lambda Orionis. The color table from Figure 1 is used, after a la 
statistical stretch on the original skyflux values. 
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